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Abstract 

A fundamental problem in plasma physics, space science, and astrophysics is the 
transport of energetic particles interacting with stochastic magnetic fields. In partic¬ 
ular the motion of particles across a large scale magnetic held is difficult to describe 
analytically. However, progress has been achieved in the recent years due to the de¬ 
velopment of the unified non-linear transport theory which can be used to describe 
magnetic held line diffusion as well as perpendicular diffusion of energetic particles. 
The latter theory agrees very well with different independently performed test- 
particle simulations. However, the theory is still based on different approximations 
and assumptions. In the current article we extend the theory by taking into account 
the hnite gyroradius of the particle motion and calculate corrections in different 
asymptotic limits. We consider different turbulence models as examples such as the 
slab model, noisy slab turbulence, and the two-dimensional model. Whereas there 
are no hnite gyroradius corrections for slab turbulence, the perpendicular diffusion 
coefficient is reduced in the other two cases. The matter investigated in this article 
is also related to the parameter ”a 1 2 ” occurring in non-linear diffusion theories. 
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1 Introduction 


Energetic particles interact with turbulent magnetic fields while they propa¬ 
gate through a magnetized plasma. Due to this interaction, their motion is a 
stochastic motion. In addition to turbulent fields 6B , there are also large scale 
magnetic fields B 0 influencing the particle orbit. The latter field is usually 
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called guide field, mean field, background field, or large scale field. This type 
of configuration can be found in different physical environments such as fusion 
devices, the solar wind, or the interstellar medium (see, e.g., Schlickeiser 2002, 
Spatschek 2008, and Shalchi 2009 for reviews). The turbulent fields described 
above lead to a diffusive motion of the particles. Due to the large scale field 
one has to distinguish between diffusion along and across that field. It is often 
assumed that perpendicular diffusion is more difficult to describe analytically 
compared to parallel transport. It should be noted that sub- and superdif- 
fusive transport has been discussed more recently in the literature (see, e.g., 
Zimbardo et al. 2006, Pommois et al. 2007, Shalchi & Kourakis 2007, and 
Zimbardo et al. 2012). 

Previous theories for perpendicular diffusion such as the Non-Linear Guiding 
Center (NLGC) theory of Matthaeus et al. (2003) and the Unified Non-Linear 
Transport (UNLT) theory of Shalchi (2010) neglect the rotation of the particle 
in the direction perpendicular with respect to the large scale field. Further¬ 
more, as equation of motion in such theories, the following Ansatz is used 

V x = av z SB x [x(t),t] /B 0 . (1) 


Here we have used the particle position x(t) at time t, the ^-component of 
the particle velocity vector v z , the x-component of the turbulent magnetic 
field 5B X , the mean magnetic field B 0 , and the x-component of the guiding 
center velocity vector V x (see Section 2 for more details). The parameter ”a” 
used here can be seen as an unknown parameter. If non-linear theories are 
compared with test-particle simulations, best agreement is usually found for 
values between a 1 2 = 1/3 and a 2 = 1 (see, e.g., Matthaeus et al. 2003 and 
Tautz & Shalchi 2011). However, analytical work based on the Newton-Lorentz 
equation suggests that this parameter is between a 2 = 1 and a 2 = 2 (see 
Shalchi & Dosch 2008, Dosch & Shalchi 2009, and Dosch et al. 2013). 

Especially the UNLT theory has shown remarkable agreement with test-particle 
simulations for different turbulence models such as the slab/2D model, the 
Goldreich-Sridhar model, Alfven waves, and noisy reduced MHD turbulence 
(see Tautz & Shalchi 2011, Shalchi 2013, Hussein & Shalchi 2014, Shalchi & 
Hussein 2014, Shalchi & Hussein 2015). Therefore, we conclude that the main 
problem in the theory of perpendicular transport is Eq. (HD and therewith the 
parameter ”a”. 

In the current paper we explore the influence of finite gyroradius effects. To 
do this we employ two different approaches, namely 

(1) Quasi-Linear Theory (QLT, see Jokipii 1966), 

(2) Non-linear diffusion theory (see, e.g., Shalchi et al. 2004, Shalchi 2010); 
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For both approaches we compute the perpendicular diffusion coefficient as a 
function of the gyroradius. Results are compared with each other and pre¬ 
vious findings are recovered by considering appropriate limits. Furthermore, 
we estimate the value of "a" in the context of finite gyroradius corrections. 
As examples we consider three different turbulence models namely the slab 
model, a noisy slab model, and the two-dimensional model. 

The remainder of this paper is organized as follows. In Section 2 we briefly 
discuss the relation between particle and guiding center coordinates and the 
corresponding diffusion coefficients. Quasi-linear theory is employed in Section 
3 and in Section 4 we use a more advanced approach based on non-linear 
diffusion theory. We end with a short summary and some conclusions in Section 
5. 


2 Relation between guiding center and particle coordinates 


Usually one is interested in the coordinates of the charged particle interacting 
with turbulence. In the following we refer to these coordinates as particle 
coordinates and we will use the symbols x and v for particle position and 
velocity, respectively. Alternatively, one can use the coordinates R defined via 
(see, e.g., Schlickeiser 2002) 


R 


x + 


c 

qB 0 


(p X e z ) 


X + — (v x e z ) 


( 2 ) 


where U = (qB 0 )/(mc'y) is the unperturbed gyrofrequency. Here we have used 
the electric charge of the particle q, the rest mass m, the speed of light c, 
and the Lorentz factor 7 . For the velocity v of the particle we use spherical 
coordinates 


v x = v\J 1 — /1 2 cos $ 
v y = v\J 1 — p? sin $ 

v z = vp, (3) 

with the particle speed v, the pitch-angle cosine p, and the gyrophase <f>. 

In the unperturbed case, where we have by definition 5B = 0, the vector R 
corresponds to the position of the guiding center. Therefore, we call R the 
guiding center coordinates. To obtain the velocity of the guiding center, we 
consider the time derivative of Eq. (J2J) . This gives us 
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( 4 ) 


-> dR ^ c f dp _ N 

V:= -dt =v + Wo[di* e ‘ 


_ 1 

= V + ~R 
-DO 


Vz B 0 V B 0 


(v x B^j 
JB Z 


x e 2 


where we have employed the Newton-Lorentz equation dp/dt = q(v x B)/c and 
the Grafimann Identity 


x (b x cj —b(a- c) — c (a -b^ . 


( 5 ) 


Therefore, we can express the components V) of the guiding center velocity 
vector by the components of the particle velocity vector Vi and the magnetic 
held components. Very often in diffusion theory, turbulence models with 5B Z = 
0 are considered. Examples are the slab and the two-dimensional model (see 
below for a definition of these two models). In this particular case Eq. (J4|) 
simplifies to 


Vr = V : 


SB x _ 

Bn 


T/ _ SBy 

Vy ~ Vz d : 

-Dn 


and V z = v z . 


( 6 ) 


From a more practical point of view, these equations are valid as long as the 
condition 5B Z B 0 is satisfied. In this case they can be used as starting point 
to compute the perpendicular diffusion coefficient. 

We have to be very careful if particle coordinates or guiding center coordinates 
are used. The guiding center coordinates (A", Y, Z ) are related to particle 
coordinates (x, y, z) via 

x = x + Y = y - and Z = z (7) 


as derived from Eq. (j2]h 

The perpendicular diffusion coefficients can be calculated as time derivatives 
of the corresponding mean square displacements. From Eq. (J7J) we find the 
following relations 


((AX) 2 ) = ((Aa;) 2 ) + ^ (AxAv y ) 
+ ^(( A ^) 2 )> 
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( 8 ) 


{(AY) 2 ) = {(Ay) 2 )-^(AyAv x ) 

+ ^((A t >,) 2 ); 

The mean square displacements ((AX) 2 ), ((AY') 2 ), ((Ax) 2 ), and ((Ay) 2 ) 
should increase linearly with time if the transport is indeed diffusive. The 
terms (AxAv y ) and (AyAv x ) correspond basically to the drift coefficient (see, 
e.g., Shalchi 2011) and, thus, they should approach asymptotically a constant. 
Since velocities are limited due to < v and \v y \ < v , the last terms in Eq. 
(|S]) can also not increase with time. Therefore, we conclude that 

((AX) 2 ) = ((Ax) 2 ) and ((AT') 2 ) = ((Ay) 2 ) (9) 

in the limit of late times. Thus the perpendicular diffusion coefficients of par¬ 
ticles and guiding centers should be the same in that limit. 


3 Quasi-linear perpendicular diffusion 


As a first step we consider Quasi-Linear Theory (QLT) and take into account 
the gyrorotation of the particle. By doing this, we basically follow Schlickeiser 
(2002). As explained in Shalchi (2015), QLT is correctly describing perpen¬ 
dicular diffusion only for very long parallel mean free paths (corresponding to 
suppressed pitch-angle scattering) and small Kubo numbers. 

For axis-symmetric turbulence, we can calculate the Fokker-Planck coefficient 
of perpendicular diffusion D± by using the Taylor-Green-Kubo formulation 
(see Taylor 1922, Green 1951, Kubo 1957) 

OO 

D± = J dt (V x (t)V x (0 )). (10) 

o 

By assuming 5B Z B 0 , we can use the first formula of Eq. OS]) to derive 

OO 

D± = lP* t J dt (v z (t)v z (0)SB x (t)5B x (0)) (11) 

0 o 

where 5B x (t) stands for SB x [x(t),t]. Within the quasi-linear approximation 
we assume that pitch-angle scattering is suppressed and, thus, 

v z (t) v z( 0) = v 2 y 2 = constant. (12) 
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Therewith, we derive 

2 2 oo 

D ± = -Jpr$l [ dt (5B x (t)5B x { 0)). (13) 

o 

For the magnetic correlation function we can employ 

(6B x (t)SB x ( 0)) = J d 3 k P xx (ky^V~ s ^). (14) 

Here we have used the magnetic correlation tensor in the Fourier space Pimik) = 
(dBi^SB^k)) and x u (t ) denotes the unperturbed orbit. Furthermore, we 
have assumed that the turbulence is static. In the following, we use cylindrical 
coordinates for the wavevector k which are related to Cartesian coordinates 
via 


k\\=k z , 

k± = yjkl + k*, 

T=arccot (k x /k y ); (15) 

According to Eq. (12.2.2a) of Schlickeiser (2002), we have within QLT 


0 ik-x u {t) _ 


+oo 


= E MW)e 


iv[ik || 


(16) 


with 

W = R L k L yJ lTF 


(17) 


and the initial gyrophase $o- The parameter Ri denotes the unperturbed 
Larrnor radius / gyroradius at fj, = 0. In Eq. (fT6|) we have also used the Bessel 
function J n (W). At the initial time t — 0, Eq. becomes 

_ +oo 

e iks u (o) = . (18) 

m =—oo 


By combining Eqs. (fill and (115]) . by averaging over the initial gyrophase, and 
by employing 

2tt 

J <Af> 0 e ^o( m -n) = 2vr(5nm , ( 19 ) 

0 
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we derive 


ifc-[£u(t)—x u ( 0 )] 


V J^(W)e i ( vllk < l+nn y. 


( 20 ) 


With Eq. (l20li we can write Eq. o as 
(5B x (t)8B x ( 0)) 

= J d 3 k P xx (k) 

x E ^(^)e < ( w/ifc ii +nn ) t (21) 

n=—oo 

and the Fokker-Planck coefficient of perpendicular diffusion (fl3|) can be writ¬ 
ten as 

D ± = n(^j 2 J d 3 kP xx (k ) 

-Too 

x E Jl{W)8(yiik\\ + nQ) 

n=—oo 

where we have used the relation (see, e.g., Hoskins 2009) 

OO 

J dt e ixt = 7 tS (x) 

0 

with Dirac’s delta distribution 5(x). In order to evaluate Eq. ()22lh one has to 
specify the tensor component P xx or consider special limits. This is done in 
the following paragraphs. 


( 22 ) 


(23) 


3.1 Zero gyroradius limit 


Eq. (j22]) depends on the unperturbed gyroradius Rl due to W = Ri,k±yJ 1 — /A 
in the Bessel functions. If we consider the limit Rl 0, and therewith W —> 0, 
we can use (see, e.g., Abramovitz & Stegun 1974) 

Jn( 0 ) = <W (24) 


and Eq. ()22li becomes 

D ± = ir J d 3 k P xx (k)5 [y/ik^j . 


(25) 
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Now we employ the relation (see, e.g., Hoskins 2009) 


5 (ax) = —5 ( x ) 
a 


(26) 


to obtain 

D ± = J d 3 k P xx (k)5 (fc||) . (27) 

In the latter formula we find the parallel integral scale Ln which is defined via 
(see, e.g., Shalchi 2014) 

5B 2 x L\\ = njd 3 k P xx (k) 5 (fc||) . (28) 


By using Ly, we can write the Fokker-Planck coefficient of perpendicular dif¬ 
fusion as 

D± = h M L, S -P-. (29) 

2 ± 5 0 

By employing (see, e.g., Schlickeiser 2002) 

1 V 

«l = 2 J D X ifi) , (30) 


we can easily calculate the perpendicular diffusion coefficient and we find 
the well-known quasi-linear result 



SB 2 

~W 


(31) 


We would like to emphasize that the latter formula is correct within the quasi- 
linear approximation and if we consider the limit Rj J —» 0. Eq. (13Tj) is usually 
called the (quasi-linear) held line random walk limit and was originally ob¬ 
tained in Jokipii (1966). 


3.2 Slab turbulence 


In the following we consider the so-called slab model for which the turbu¬ 
lent magnetic held depends only on the coordinate along the guide held, i.e., 





8B{x) = SB(z). In this particular case the components of the magnetic corre¬ 
lation tensor are given by 


Pim(k) 




Im • 


(32) 


Here we have used Kronecker’s delta 5i m , Dirac’s delta distribution 8(z ), and 
l,m = x, y. To satisfy the solenoidal constraint, we need to have Pi z = P zm = 
P zz = 0. Here g slab {k\\) is the turbulence spectrum of the slab modes. 

If we combine Eqs. (I22|) and (1321) we find 
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D± = 2tt 2 J dk, g M {k^ 

—oo 

+oo 

x X (yp^w + n • 


Employing again Eq. (1241) is leading to 

Di = 2 jt 2 (|h J dk, g sLab (k,)6 (fMn) , 

—OO 

Now we use Eq. to obtain 


(33) 


(34) 



—OO 


g‘ m (k,)S (k,) 


r, 2^ 1/^1 „slab(i 
= 2 

To continue, we need to 
following we employ 

g slab {k\\) = ^ySB 2 l\ 


= 0 ). 

specify the function g slab (k \\) 
r , . oi —s/2 

i + (fc||/n) 


(35) 

in Eq. (I35|h In the 

(36) 


as suggested by Bieber et al. (1994). Here we have used the strength of the 
turbulent held SB and the bendover scale l\\ indicating the turnover from the 
energy range of the spectrum to the inertial range. The normalization function 
C(s ) = T(s/2)/[2i/Tr((s —1)/2)] depends on the inertial range spectral index 
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s and contains Gamma functions. Therewith, Eq. (j35j) becomes 


D ± = ttC(s)v \h\ l\\ —^2 


'o 


( 37 ) 


One can easily show that for the spectrum used here, the two parallel length 
scales are related to each other via Ly = 27rC'(s)Zy. Therefore, Eq. (T37l) agrees 
with Eq. 029]) and, thus, Eq. (13T1) is also valid in the case considered here. 

We have shown that for slab turbulence, the zero gyroradius limit is exact if 
quasi-linear theory is employed. This result was predictable because for slab 
turbulence we have by definition 5B X = 5B x (z ) and, thus, the perpendicular 
motion of the particle is not relevant. 

3.3 A noisy slab model 

The slab model used in the previous paragraph is an extreme model due to 
the fact that there is absolutely no transverse structure in the turbulence. 
Therefore, we extend the slab model by employing 


2 / 

Pm, (k) = jVyfcnje (i - kj ± ) 



(38) 


where we have used the Heaviside step function Q(x). We refer to this model 
as the noisy slab model which was originally described in Shalchi (2015). Due 
to the step function, there is only turbulence if k±l± < 1. For large enough 
wavenumbers, on the other hand, there is no turbulence. Physically this cor¬ 
responds to the case where wave vectors are mainly oriented parallel with 
respect to the mean field but weak fluctuations are taken into account. The 
idea of incorporating some kind of noisiness via a step function goes back 
to Ruffolo & Matthaeus (2013) where this idea was combined with the two- 
dimensional turbulence model. A more detailed explanation of this type of 
turbulence model can be found there. 

Combining the noisy slab model represented by Eq. (138|) with Eq. (1221) yields 



-1 
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X 


(39) 


—(—oo 

J dk\\ g slab {k\\)5(ynk\\+nti) . 

— OO 


Now we employ Eq. (1261) and assume that the spectrum is symmetric, i.e., 
g slab {k\\) = g slab (—k\\). After evaluating the fey-integral we obtain 


D I = 27 T l 


2, V \^\ 


±- 


Bl 


dk 


x 


+oo 


^ E Jn(w)g 


slab 



To proceed we employ spectrum (136|) to deduce 


(40) 


D±_ = irC(s)/||/j.^-t) \p\ 


+oo 


X 


dk _l J2 


J 2 n(W) 


[ n 2 + (fiRY 


s/2 


(41) 


where we have used R = Rl/ 1\\ = v /{^\\)- We can easily recover the pure slab 
result (137)1 by considering the limit l± —> oo. 

To achieve a further simplification of Eq. (ITTlh we consider the special case 
s = 2. In this case we can use C(s = 2 ) = 1/(27t) to write 


D± = h\\l±^v\g,\ (/iR) 2 


+oo 


x / dk± 


J 2 n(W) 


=-00 n 2 + (fiR) 


2 • 


(42) 


We would like to add, that in this case the parallel integral scale Ly and the 
bendover scale Zy are equal. In order to evaluate the form (I42p analytically, 
one has to approximate the series involving Bessel functions. The latter series 
is discussed in detail in the Appendix of the current paper. These discussions 
are based on the work of Shalchi & Schlickeiser (2004) and Tautz & Lerche 
(2010), respectively. 
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3-4 Noisy slab turbulence and small gyroradii 


In the following we focus on the case of small gyroradius corrections. Therefore, 
we assume that the parameters W = Ri,k±^l — /x 2 and fiR are small as well. 
In this particular case we can employ Eq. (I94[) . Therewith, Eq. (j42|) can be 
approximated by 


D ] 


h, S J1 

2 1 Bi 


■v Ihl 


r r l 

X j dl:_ 1 - - (l - y) K 2 L k j 


(43) 


The remaining wavenumber integrals can easily be solved and we deduce 
1 SB 2 


D i 


i l " Bl 


■V \fx\ 


l r >\ R\ 




l 2 


(44) 


Employing Eq. (1301) allows as to compute the perpendicular diffusion coeffi¬ 
cient. After straightforward algebra we fold 





SB 2 

AT 




(45) 


where we have used again R = Rl/1\\- Or, if one is interested in the perpen¬ 
dicular mean free path 


Ax 


3 SB 2 

4 "TT 




(46) 


If the parameter ”a 2 ” used in Eq. (JT[) is understood as finite Larmor radius 
corrections, we ford for the case considered here 


1 ll 


1--f R z . 

12l\ 


(47) 


Here we can directly see, that finite gyroradius corrections depend on the nor¬ 
malized gyroradius R and the scale ratio l\\/l±- This conclusion is in agreement 
with Fig. [0 Furthermore, we conclude that the parameter ”a 2 ” is smaller than 
one. 
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Fig. 1. The perpendicular mean free path Aj_ versus the unperturbed gyroradius 
Rl■ Both parameters are normalized with respect to the parallel scale Zy. All curves 
shown here are obtained for quasi-linear transport in noisy slab turbulence. The 
solid line corresponds to the zero Larmor radius limit and the other curves are 
obtained by taking into account the non-vanishing gyroradius. Shown are the results 
for Z||/Z_l = 0.1 (dotted line), l\\/l± = 1 (dashed line), and = 10 (dash-dotted 
line). 

3.5 Numerical results for noisy slab turbulence 



For arbitrary Larmor radii Rl , we can compute the perpendicular mean free 
path only by solving Eq. (lT2j) numerically. The results are visualized in Fig. [0 
There, we have also shown the zero Larmor-radius limit. The latter limit can 
be obtained from Eq. (1461) by setting R = 0 therein. In this particular case 
the perpendicular mean free path becomes 


A_l _ 3 SB 2 


(48) 


Furthermore, we have computed the ratio X±/l\\ versus Rl/1\\ f° r different 
values of the ratio l\\/l±■ In all considered cases we find that finite Larmor- 
radius corrections reduce the perpendicular diffusion coefficient. For l± = 0.1/||, 
for instance, the perpendicular mean free path is reduced by a factor 3. 


3.6 Two-dimensional turbulence 


Another popular model for turbulence is the two-dimensional (2D) model (see, 
e.g., Shalchi 2009 for details). In this case the magnetic correlation tensor is 
given by 

Pt,a(k) = g 2D (k ± ) S -lRi (* m - ^ (49) 


13 









if l , m = x,y. The tensor used here corresponds to the one dehned in Bieber 
et al. (1994). In the latter paper it was also assumed that 8B Z = 0. Therefore, 
tensor components with l = z and/or m = z are zero. The function g 2D (k±) 
used above is the spectrum of the two-dimensional modes. 

By combining Eq. (149|) with Eq. (1221) we find after straightforward algebra 


D j_ = 7r 2 



J dk± g 2D {k±) 

o 


x 


+oo 

E WM- 


n=—oo 


(50) 


Only the contribution with n = 0 is not zero due to the Dirac delta 5(nQ). 
For n — 0, however, we have 5(0) = oo and, thus, D± = oo. Therefore, QLT 
provides for the two-dimensional turbulence model an infinite perpendicular 
diffusion coefficient. This result is unphysical and was found before in diffusion 
theory (see, e.g., Shalchi et al. 2004). Therefore, QLT cannot be used in order 
to describe perpendicular transport in the general case. In the next section we 
consider transport in two-dimensional turbulence again but in the context of 
non-linear diffusion theory. 


4 Non-linear diffusion theory 


It is well-known and accepted that quasi-linear theory cannot describe per¬ 
pendicular diffusion in the general case (see, e.g., Spatschek 2008 and Shalchi 
2009). Even if one would assume that pitch-angle scattering is suppressed and 
perpendicular diffusion is only caused due to the stochasticity of magnetic field 
lines, quasi-linear theory would not work due to the fact that the field line ran¬ 
dom walk itself is non-linear for large Kubo numbers (see, e.g., Shalchi 2015). 
Therefore, a non-linear treatment is required in order to describe perpendicular 
transport. A very powerful tool is the Uuihed-Non-Linear Transport (UNLT) 
theory of Shalchi (2010). The latter theory can describe transport of magnetic 
field lines and particles across the mean field. The theories of Matthaeus et al. 
(1995) and Matthaeus et al. (2003) as well as quasi-linear theory are contained 
in this theory as special limits. 

In the following we use ideas of non-linear diffusion theory for the special case 
of suppressed pitch-angle scattering. The case of strong pitch-angle scattering 
will be investigated in the sequel of the current paper. According to Eq. dill) 
the critical quantity in the theory of perpendicular diffusion is the 4th order 


14 


correlation function 


C{t 1 ,t 2 ) := ( v z (ti)v z (t 2 )SB x (ti)8B x (t 2 )) 


(51) 


where we have used again the notation 5B x (t) = ARja^f), t]. Such 4th order 
correlations are difficult to evaluate. Matthaeus et al. (2003), for instance, 
suggested to approximate 4th order correlations by a product of two 2nd order 
correlations. In the general case, however, this approximation is not valid (see 
Shalchi 2010). In the current paper we consider only the case that velocity 
diffusion is suppressed and, thus, we can approximate 


{v z (t)v z (0)5B x (t)5B x (0)) 

« uV J d 3 k P xx (k) . (52) 

Here we have also employed the so-called random-phase approximation which 
is a well-known tool in diffusion theory (see, e.g., Lerche 1973, Matthaeus et al. 
1995, Matthaeus et al. 2003, Tautz & Shalchi 2010). The latter approximation 
is sometimes called Corrsin’s independence hypothesis (see Corrsin 1959). 

Here one has to be careful since the magnetic held has to be computed at 
the particle position. Therefore, x{t) and (7(0) are particle positions and not 
guiding center positions. To evaluate the characteristic function (exp (ik ■ x)), 
previous authors used the propagator of the diffusion equation (see Matthaeus 
et al. 2003) or an approach based on the Fokker-Planck equation (see Shalchi 
2010). However, the latter two equations are derived by using guiding center 
coordinates (see, e.g., Schlickeiser 2002). 

With Eq. (17jl . the characteristic function in Eq. (l52]i becomes 


Jk-[x(t)-x( o)]^ = / ik-[x(t)-X( 0)]\ 


ikxAvy/Q 

Alternatively, one can write 


^ik-[x(t)—x{ 0)] \ __ ^ik-[x u (t)— x u (0)] PiAXk x +iAYk y 


(53) 


(54) 


where we have used the unperturbed orbit x u (t ) as in Section 3. 

The UNLT theory allows in principle to calculate D± non-lincarly. By using the 
latter approach, 4th order correlations of the form (v x (t)v x (0)5B x (t)5B x (0)) 
can be computed by using the Fokker-Planck equation which has the form 
(see, e.g., Schlickeiser 2002) 
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( 55 ) 


df , df df 

— + V/JL — ~ \l — 

^ dZ d <P 





Here we have assumed axis-symmetry with respect to the z-axis and we have 
used the Fokker-Planck coefficients of pitch-angle diffusion D w and gyrophase 
diffusion respectively. In non-linear diffusion theories a fundamental 

quantity is the so-called characteristic function which is defined as 

r (k, p, t) := f d 3 X e*-*f(X, /i, II 0 , t) = (e*-*) . (56) 


To proceed we multiply the Fokker-Planck equation ([55]) by exp (ik ■ A") and 
thereafter we integrate over space to obtain 


ar 

~dt 

d_ 

d/i 


df 


“a* 


r_ an 

d 

' ar' 

i- 

b 

■C 

-§1 

+ a¥ 

a$ 


D ± k 2 ± T. 


(57) 


If we suppress pitch-angle scattering and gyrophase diffusion by setting D /t/t = 
0 and Dq,q, = 0, the latter equation is solved by 

r (k,(jL,t) = e ikzVIMt - D±k i t (58) 


where we have also average over the gyrophase <f>. This formula provides the 
characteristic function for suppressed velocity diffusion. Therefore, Eq. (15TD 
becomes in the case considered here 


^ik-[x(t)-x{ 0 )]^ _ e ik-[x u (t)—x u (0)] e —D±k\t 


(59) 


where x u (t) denotes the unperturbed orbit as before. 

To proceed we combine Eqs. and (l59j) to derive for the Fokker- 

Planck coefficient of perpendicular diffusion 

v 2 if 2 +°° r - 

D2 - = -XT E d 3 k R n (k)P xx (k)4(W) (60) 

-°0 n=—oo J 

with the resonance function 
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R n (k) — Re J dt e i{v ^ +nn)t e- D ^ 

o 

D±kj 

(D±k'Y) 2 + {vjik\\ + nil) 


( 61 ) 


We can easily see that perpendicular diffusion described by the parameter D± 
broadens the resonance. The effect of resonance broadening due to perpendic¬ 
ular diffusion was already described in Shalchi et al. (2004). 


The quasi-linear limit (122]) can be recovered by considering D± 0 and by 
using (see, e.g., Hoskins 2009) 


lim —-- 

a->o a -f- x- 


7 t5 ( X ) . 


(62) 


This is leading to the results discussed in Section 3 of the present paper. 

In the current section, Eq. (1661) with (16TT) was derived by using the approach 
proposed in Shalchi (2010) for a finite gyroradius but suppressed velocity 
diffusion. We would like to emphasize that the same result can be obtained 
from the so-called Weakly Non-Linear Theory (WNLT) of Shalchi et al. (2004) 
if velocity diffusion is suppressed therein. 


4-1 Zero gyroradius limit 


Eq. (1661) with (T6TT) describes perpendicular transport for suppressed velocity 
diffusion but with finite gyorradius in the non-linear case. We can recover the 
zero gyroradius limit by considering W —> 0. In this case we can use again Eq. 
(1241) and Eq. (l60l) becomes 

2 2 

Dx= l ^r j d 3 k R 0 (k)P xx (k ) (63) 

with 

uk = ——v- < 64 > 

{D ± kD + (vfj,k\\j 

This result agrees with the UNLT theory for the case of suppressed pitch-angle 
scattering. As already discussed in Shalchi (2010), Eq. (1651) with (1641) can be 
solved by 

D±((jl) = v \n\ k F l (65) 
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where Kfl is the solution of the following non-linear integral equation 


Hfl — 




drk P x 


(k) 


KFLk'] 


(liFLkl) + kl 


( 66 ) 


Eq. (165|) corresponds to the Field Line Random Walk (FLRW) limit and the 
parameter kfl is the diffusion coefficient of field line wandering. Eq. (I66p is a 
non-linear integral equation for kfl which was originally derived by Matthaeus 
et al. (1995). 


4-2 Slab turbulence 


In the following we employ again the slab model defined via Eq. (j32]h If the 
latter model is combined with Eq. (I60jh we deduce 


?; 2,,2 +oo 

d±= 2* v -±- y m 

-^0 n =—oo 

+00 


x J dk\\ g sla \k l] )R n (k lll k ± = 0). 


To proceed we use again Eq. (1241) to obtain 
+ 00 


D i =2 tt 


U 2 /i 2 


dk\\ g slab (k ll )R 0 (k h k ± = 0). 


(67) 


( 68 ) 


Now we combine Eq. (164|) with (162]) . and we employ again Eq. (126|) to fold 

I I + 00 


D± = / dk {{ g slab (k\\)5(k\\) 

n 0 J 


(69) 


which agrees with the slab result derived above. Obviously, finite Larmor ra¬ 
dius corrections are not important in this case. This result was predictable, 
because in the current section we only included the factor T(fc|| = 0) = 
exp(— D±k‘j_t) to the quasi-linear formula. For slab turbulence this factor is 
one. 


18 






4-3 Two-dimensional turbulence 


As shown above, the non-linear effect due to perpendicular diffusion is not 
relevant for slab turbulence. Therefore, we have to consider a model with 
transverse structure. A simple model, for which non-linear effects are impor¬ 
tant, is the two-dimensional model used already above in the context of QLT. 
For this turbulence model we can employ Eq. (1491) and therewith Eq. (1601) 
becomes 


2 2 +oo °? 

D± = n J2 f dk ±g 2D (k ± ) 

-°0 n=—oo J 0 


x R n (fc|, = 0, k ± ) J 2 n (W) 

(70) 

with 


R n (k\\ = 0, k±) = -—-o • 

V J ( D ± k 2 ± f + (nil) 2 

(71) 


As shown in Shalchi et al. (2004), the latter formulas for suppressed velocity 
diffusion can be written as 

2 2 oo 

Dl = / dk ± g 2D (k ± )VH(V, W) (72) 

0 0 


where we have used the series 


H(V, W) 


+oo 

E 


n=—oo 


V 2 + n 2 


(73) 


and the parameter V = D±k'j_/Q. The series ()73j) occurred already above in 
the context of QLT (see Eq. (142)1 ). A more detailed discussion of this series 
can be found in the Appendix of the current article. 

For the spectrum of the two-dimensional modes, we employ the model pro¬ 
posed by Shalchi & Weinhorst (2009) 


~ 2D “- ' 2D ^ q) -8B% 


g (k±) = 


7T 


X 


(k±l±Y 


[1 + (k ± l ± ) 2 ] {s+q)/2 ' 


(74) 
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Here we have used the perpendicular bendover scale l±, the inertial range 
spectral index s, and the energy range spectral index q. For the inertial range 
spectral index we employ s = 5/3 as suggested by Kolmogorov (1941). The 
physical consequences of the different values of q are discussed in Matthaeus et 
al. (2007) and Shalchi & Weinhorst (2009). It was shown there that in order 
to obtain a finite ultra-scale, the condition q > 1 must be satisfied. In the 
current paper we set q = 3 as an example. In the spectrum (1741) we have used 
the normalization function 

2r (i=!) r (s±i) 



with the Gamma function T(z). 


4-4 Two-dimensional turbulence and small gyroradii 


In the following we evaluate Eq. (1721) in the limit of small gyroradii. We do 
this by considering the asymptotic limit Rl —>■ 0. Therefore, the parameter 
W, which is directly proportional to Rl, can be assumed to be small. More 
complicated to estimate is the parameter V which can be written as 


V 


D±kj K±k 2 \ 

ft ~ ft 


1 Rl A_l 
3 l_ i_ 11 _ 


(Lk ± ) 2 . 


(76) 


Since we consider the limit Rl —> 0, we can also assume that V is small, at least 
as long as the perpendicular mean free path is smaller than the perpendicular 
scale l±. The latter condition is usually satisfied (see, e.g., Fig. [2] of the current 
paper). Thus, we have to approximate the series (1731) for the case of small 
values of V and W. In this special case we can employ approximation (194)) 
and Eq. (1771) becomes 


D \ ~ 7T 



J dk± g 2D (k±)k ± 2 

o 


X 


1 



(l - II 2 ) k\ 


(77) 


The first integral can be replaced by the diffusion coefficient of random walking 
magnetic field lines (see, e.g., Shalchi 2014) 


k fl = J 2 / dk± 9 2D (k±)k ± 2 
0 0 


(78) 


20 







and the second one by the normalization condition (see, e.g., Shalchi 2009) 


CXJ 

J dk ± g 2D (k_ l) = 


SB 2 

~27 


(79) 


Therewith, Eq. (1771) becomes 


D\ 


!‘i>f k'fl - \ m 2 (i - y) R i-§T' 


(80) 


We can easily take the square root of the latter equation. By assuming that the 
second term in Eq. (I80p is much smaller than the Erst one, we can approximate 
the Fokker Planck coefficient of perpendicular diffusion by 


D, 


v |/i| k fl 


, 1 o' 

R\ 

SB 2 1 


K 2 

K FL 

Bl\ 


(81) 


The first term corresponds to the well-known held line random walk limit. 
In order to compute the perpendicular diffusion coefficient k±, we have to 
combine Eqs. ()30l) and (j8T]l to End 


v 

k± ~ -kfl 


1 R\ SB 2 

16 Hp L Bq 


(82) 


For spectrum (1741) the field line difiusion coefficient was calculated in Shalchi 
& Weinhorst (2009) based on the Matthaeus et al. (1995) theory. They found 


s — 1 5B 


(83) 


as long as q > 1 is satisfied. By combining Eq. (1831) with (1821) . we can easily 
compute the perpendicular diffusion coefficient or the perpendicular mean free 
path A_l = 3k±/v. We derive 


Au 


3 

2 Kfl 



q- i Rl\ 
8( S - l) (i j' 


(84) 


In this particular case the parameter ”a 2 ” occurring in non-linear diffusion 
theory has the value 


g~l R 2 l 

8(s — 1) l\ ' 


(85) 
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Fig. 2. The perpendicular mean free path Aj_ versus the unperturbed gyroradius Rl- 
The first parameter is normalized with respect to the perpendicular scale l± and 
the second parameter to the parallel scale Zy. All curves shown here are obtained 
for non-linear transport in two-dimensional turbulence. We set the energy range 
spectral index q = 3 to ensure that all turbulence scales are finite. The solid line 
corresponds to the zero Larrnor radius limit and the other curves are obtained 
by taking into account the non-vanishing gyroradius. Shown are the results for 
Zi|/Z_l = 0.1 (dotted line), Zm/Zj_ = 1 (dashed line), and Z||/Zj_ = 10 (dash-dotted 
line). 

Again we found that finite gyroradius effects reduce the perpendicular diffusion 
coefficient and, thus, a 2 < 1. Furthermore, the parameter ”a 2 ” depends on the 
scale l± and the two spectral indexes s and q. 


4-5 Numerical results for two-dimensional turbulence 


For arbitrary Larrnor radii Rl, we can compute the perpendicular mean free 
path only by solving Eq. (1721) numerically. The results are visualized in Fig. 
[2] There, we have also shown the zero Larmor-radius result. The latter limit 
can be obtained from Eq. (1841) by setting Rl = 0 therein. In this case the 
perpendicular mean free path becomes 

Aj_ _ 3 hfl _ 3 / s — 1 SB 

77 “ 217 ~~ 2 y 2(q — 1) 77 1 j 


where we have employed Eq. (1831) also. For s = 5/3, q = 3, and SB = B 0 we 
find A_i_/Z_i_ ~ 0.612. Furthermore, we have computed Aj_/Zj_ versus Rl/1\\ for 
different values of the ratio Z||/Zj_. In all considered cases we find that finite 
Larmor-radius corrections reduce the perpendicular diffusion coefficient. For 
lj_ = 0.1Z||, for instance, the perpendicular mean free path is about a factor 
3 smaller compared to the zero Larmor-radius limit. If the parameter ”a 2 ” is 
understood as finite gyroradius effect, this means that a 2 ~ 1/3 for the case 
l _l = 0.1Z||. That is exactly in agreement with what was obtained by Matthaeus 
et al. (2003). 
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5 Summary and conclusion 


In the current paper we have investigated the effect of a finite gyroradius in 
the analytical theory of perpendicular diffusion. To achieve a simplification of 
our calculations, we suppressed the non-linear effect of velocity diffusion. The 
case of strong velocity diffusion will be explored in the second paper of the 
series. 

In order to better understand finite gyroradius effects, we have combined two 
analytical theories for perpendicular diffusion with different turbulence mod¬ 
els. We considered quasi-lincar and non-linear transport, respectively. For the 
turbulence we have employed the slab model, a noisy slab model, and the 
two-dimensional model. The first model represents the case that the Kubo 
number is zero, whereas the other two models represent turbulence with small 
and infinite Kubo numbers, respectively. 

We have shown that finite gyroradius effects do not occur in slab turbulence as 
expected. For the other two models such effects can be important depending 
on the size of the gyroradius. In all considered cases, finite gyroradius effects 
reduce the perpendicular diffusion coefficient. A reduction of the perpendicular 
diffusion coefficient due to finite gyroradius effects was also found in Neuer & 
Spatschek (2006). 

A reduction of the perpendicular diffusion coefficient can be important for 
describing particle acceleration at shock waves (see, e.g., Ferrand et al. 2014). 
In such scenarios a small diffusion coefficient leads to a longer residence time of 
the particles at the shock front and, thus, the particles experience acceleration 
to higher energies. 

A fundamental problem in analytical treatments of perpendicular transport 
is the physical meaning and the value of the parameter ”a” as used in Eq. 
©. One explanation of this parameter is provided in the current paper. The 
parameter discussed here can be understood as finite gyroradius effect. In all 
cases found in the current paper, the finite gyroradius reduces the perpendicu¬ 
lar diffusion coefficient corresponding to a value a 2 < 1. This is basically what 
one finds if analytical results are compared with test-particle simulations (see, 
e.g., Matthaeus et al. 2003 and Tautz & Shalchi 2011). We have also shown 
that the absolute value of ”a” depends on the turbulence length scales l\\ and 
l l as well as on the inertial range spectral index s and the energy range spec¬ 
tral index q. For two-dimensional turbulence and l± = 0.1/||, for instance, we 
found a 2 ~ 1/3. This is exactly what was found before in Matthaeus et al. 
(2003). In the sequel of the current paper we shall explore the case of strong 
velocity diffusion. 


23 


Acknowledgements 


Support by the Natural Sciences and Engineering Research Council (NSERC) 
of Canada is acknowledged. 


Appendix: The Series H(x,y) 


In the current paper the series 



(87) 


occurs in the quasi-linear treatment of the transport in noisy slab turbulence 
(see Eq. (142 h ) and in the non-linear treatment of the transport in pure two- 
dimensional turbulence (see Eq. (1731) ). Eq. (1571) is a special type of a so-called 
Kapteyn series (see Kapteyn 1893 and Watson 1966). This special type of se¬ 
ries was already discussed in Shalchi & Schlickeiser (2004) where the following 
asymptotic limits have been derived 


T 2 (x) 1 

if(z»l,2/«l)»-^l» — 



y 2 



( 88 ) 


y 2 


A very detailed and useful investigation of series (157)) can be found in Tautz 
& Lerche (2010). These authors have shown that the series can be written as 


y\ _ n J+iy( x )J-iy( x ) 


y sinh(7n/) 


(89) 


In the following we consider small and large arguments x of the Bessel func¬ 
tions and simplify Eq. (1591) in these cases. 
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The limit Id —> 0 


In the limit considered here, the Bessel functions can be approximated by a 
Taylor expansion (see, e.g., Watson 1966) 


oo (-l) n 

Ju ^ S n \r {v + n + 1) 


\ 2n+v 

2 , 


(90) 


where we have used the Gamma function T(z). The latter expansion is con¬ 
vergent for |x| < oo and arbitrary u. Therewith, we can show that for small 
x 


J--\-iy (^) J—iy (*^) 


T(l + iy)r(l-iy) 


x 


x 


2 1 + y 2 


(91) 


where we have also used the relation T(z + 1) = zT{z). Furthermore, we can 
employ (see, e.g., Abramowitz & Stegun 1974) 


T(1 + iy)T{ 1 - iy)=iyT{iy)T(l - iy) 

—i sin(z7n/) 

ny 

sinh(7 ry) 


(92) 


With that result and by employing Eq. (l9Tj) . we can write Eq. (I89|) as 


H(x,y) « 



1 x 2 \ 
21 + y 2 )' 


(93) 


If we additionally assume that y •C 1, we can approximate 

H(x,y) « ^ (i - ^x 2 ) . (94) 

If the term which is quadratic in x is neglected, we obtain one of the limits 
listed in Eq. (155]) . If y 1, Eq. ilTIl) would become 

M « 1 (1 - W) (95) 


and in lowest order x/y we would still find the same limit as in Eq. (l88lh 
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The limit Id —> oo 


If the argument x in the Bessel function is large, we can approximate (see, 
e.g., Watson 1966) 


Jv (-^) 



(96) 


Therefore, we find 


J-\-iy (*^0 J—iy (*^) 


7T 


— cos x - ly - 

TTX 


7r 7T 

x cos ( x + —iy ~ 


By using the relation 

cos (a ± h) — cos (a) cos ( b ) =|= sin (a) sin (6) 


(97) 


(98) 


this becomes 


J-\-iy (*^) J—iy (x) 



■ 2 ( ' K \ ■ 2 ( n ■ \ 
— sm hr — — I sm I —iyJ 


which can be written as 


(99) 


J+iy (x) ^—iy (^0 


2 

TTX L 


7T 

cos“ [ —iy 


Sill 


7T 

x _ 4 


( 100 ) 


The cosine function can be replaced by using cos(iz) = cosh(^). Combining 
our findings with Eq. (l89lh we obtain 


H(x,y) w 


2 

xy sinh(7n/) 


x 



( 101 ) 
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The latter approximation can be used in the limit |x| —> oo. If we additionally 
assume that \y\ —> 0, Eq. (11011) becomes 



( 102 ) 


in agreement with Eq. (|88j) . For \y\ —> oo, however, Eq. (I101D becomes 


1 


H(x,y) 


(103) 


xy 


which is also listed in Eq. (15511 . 
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